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Abstract 

We discuss the electronic properties of graphene and graphene nanoribbons including 
"pseudo-Rashba" spin-orbit coupling. After summarizing the bulk properties, we first 
analyze the scattering behavior close to an infinite mass and zigzag boundary. For low 
energies, we observe strong deviations from the usual spin-conserving behavior at high 
energies such as reflection acting as spin polarizer or switch. This results in a spin polar- 
ization along the direction of the boundary due to the appearance of evanescent modes in 
the case of non-equilibrium or when there is no coherence between the two one-particle 
branches. We then discuss the spin and density distribution of graphene nanoribbons. 

Pacs: 71.70.Ej,73.61.Wp,72.25.Rb 

1 Introduction 

Graphene, the single-layer allotrope of carbon, is undoubtedly one of the most active 
fields in today's both experimental and theoretical condensed matter physics [H El Ej. 
Among an entire plethora of phenomena and proposals, the issue of spin-orbit coupling 
has generated particular interest [H O IH [71 |H1 • A detailed understanding of spin-orbit 
interaction in graphene is crucial for the interpretation of ongoing experiments on spin 
transport performed by various groups [TOl [TTl [121 [131 (HI [El Ell [13 [H] ■ Other issues 
include various device proposals [THl HO] and theoretical predictions dH [551 US] related to 
spins and spin-orbit coupling in graphene. 

In the present paper we investigate a single layer of graphene in the presence of spin- 
orbit interaction of the "pseudo-Rashba" type, coupling the sublattice or pseudo spin to 
the physical electron spin [H [5l [6l [7l [H [9l [24] . Our interest is based on the fact that for 
graphene on Ni with intercalation of Au a 100-fold enhancement of the "pseudo-Rashba" 
spin-orbit coupling has been reported [25j . Furthermore, impurities which induce a sp3- 
distortion will lead to a "pseudo-Rashba" spin-orbit coupling with a value comparable to 
the one found in diamond and other zinc-blende semiconductors |26j . The latter result 
indicates that the "pseudo-Rashba" spin-orbit coupling can be controlled via the impurity 
coverage. 
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In this paper, we will concentrate on the scattering behavior of spin densities near 
boundaries created either by an infinite mass or a zigzag edge. Our presentation is orga- 
nized as follows: In section [5] we introduce the basic Hamiltonian and discuss its general 
bulk solution in the absence of a mass term; the technically more complicated case of a 
nonzero mass is deferred to the appendices. In the following section |3] we investigate in 
detail the scattering properties and spin dephasing at hard boundaries for various types 
of incoming spinors and energy ranges. This discussion is extended in section [4] to aver- 
aged spin polarizations obtained from continuous distributions of incoming directions. In 
section [SJ we analyze the spin and density distribution of graphene nanoribbons. We close 
with a summary in section [6] Throughout this manuscript, we use parameters of Ref. (25] . 



2 Dirac fermions with "pseudo-Rashba" spin-orbit cou- 
pling 

The single-particle Hamiltonian of monolayer graphene with "pseudo-Rashba" spin-orbit 
interaction can be formulated as [U El O El] 



Ti, = vpp ■ T + \ {f X a) ■ ez ■ 



(1) 



where, among standard notation, A is the spin-orbit coupling parameter, and the Pauli 
matrices f , a describe the sublattice and the electron spin degree of freedom, respectively. 
For a given wave vector k this Hamiltonian reads explicitly: 
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From experience with the "classic" Dirac equation of relativistic quantum mechanics, it is 
occasionally of use not to study just a given Hamiltonian but also its square. Here we find 



H^(k) = 



{hvpk)^ 
2i\hvp{kx + iky 





-2iXhvp{kx — iky) 

{hvpk)^ + 4X^ 

invpk)^ + AX^ 

2iXhvp{kx + iky) 



This matrix is block-diagonal with eigenvalues 



{e^)i,2 = (nvpk)^ + 2X^ ± 2\X\y/ {hvpk)'^ + X'^ 
where the positive sign corresponds to the eigenvectors 



\ 



-2iXhvp{kx — iky) 
{hvpk)^ J 

(3) 



(4) 



/ sin(i9/2) \ 
cos(t9/2)e"' 








cos(i9/2 
y sin(i9/2)e*'' / 



(5) 



while for the negative sign we have 



m) = 



cos(t9/2) \ 
-sin(i9/2)e"' 





1/32) = 



-sin(?9/2) 
\ cos(?9/2)e*'' J 



(6) 
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where ^ e [0, tt] and 

C0S1? 



|A| 



A i{kx + iky) 



^{nvpkY + A2 ' |A| 
In the basis (|q;i), |a2), |/32)) the Hamiltonian reads 



k 



(7) 



/ \ 

g+ 

g_ 

\ y 



g± = ±TiVF{kx ±iky)f±{\\\/TiVFk) 



(8) 
(9) 

f±{x) = Vl + a;2±a;. (10) 
Now it is straightforward to obtain the full cigcnsystcm: Wc find a gaped pair of eigenvalues 

si,± = ±(^/{hvFkY + X' + \\\) (11) 

with eigenspinors (type I) 



with 
and 



\x.Ak)) = ^ 



sin(i?/2) \ 



cos(??/2)e''' 
±cos(i?/2)e^'^ 
\ ± sin(^?/2)e'''e*'^ / 



and 



iV _ + iky 

~ k 



With gv = 2 being the valley degeneracy, the corresponding density of states reads 



gv 



{\e\-\X\)9{e'-{2Xr) 



The other pair of dispersion branches does not exhibit a gap, 



£2,± = ± 
and has eigenspinors (type II) 



^JihvFkf + A2 - |A| 



(12) 



(13) 



(14) 



(15) 



IX2,±(fc)) 



1 



V2 



/ cos(i?/2) \ 
-sin(t?/2)e^'' 
± sin(t?/2)e*'^ 
\ T cos(t?/2)e'''e''^ / 



(16) 



The corresponding density of states reads 

P2(e) 



gv 



(N + |A|) 



2'K{hVFf 

Let us now consider expectation values within the eigenstates with wave functions 



(r1fc,/i,±) = — |x;,,±) 



(17) 



(18) 
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Figure 1: A plane wave of type I with spin perpendicular to the momentum k = {kx,ky) 
(if = avctan{ky / kx) , (po- = (p + 7r/2) is reflected at the boundary into a plane wave with 
k' = ikx,—ky) and k" = {kx,—ky) with perpendicular spin, but anti-parallel with respect to 
each other (see Eq. (j24p for the definition of ky ). 



/X e {1,2}, and A being the area of the system. Here we find 



A 



sin 1? cos if 



(A:,1,±|t|/c,1,±) = (fc,2,±|f|fc,2,±) =±— I siniJsin^ |, (19) 



■|A| 

and 



A 



■ sin 1? sin 



(A:,1,±|ct|/c,1,±) = -(A:,2,±|ct|/c,2,±) = — ( sin?9cos(p |. (20) 



|A| 

Here, if is the usual azimuthal angle of the wave vector, k — k{cosip, sin ip). Note that 

(f).(a>-fc.(a)=0, (21) 

as usual for Rashba spin-orbit coupling, and 

|(r)| = |(a)Hsin^, (22) 

where for sini? < 1 sublattice and electron spin degree of freedom are entangled which 
each other. 



3 Spin dephasing due to reflection on a hard wall 



In this section, we will study the scattering behavior from a hard wall which will lead to 
spin dephasing as depicted in Fig. [TJ For that, a general plane wave with fixed momentum 
kx and energy E > 2|A| is written as 



tpE,k, {x,y) 



with 



A,e<y\xiAk.-kl)) ^ A2e''^ly\X2Ak..kl)) 

{hvpk^f = {E + i-mMf - A' - {hvFk,)' 



(23) 
(24) 
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H G {1,2} and the normalization constant Afj:. For energies E < 2|A|, some modification 
in Eq. have to be made which shaU be discussed in more detail below. 

In the following, we will discuss the reflection at a hard wall at y = for the two types 
of plane waves, i.e., we will first set Ai = 1, A2 = (type I) and then Ai = 0, A2 — 1 
(type II). The discussion is based on the reflected spin direction which shall be denoted by 
if'^. It is obtained from the expectation value of the spin-density operator at the boundary 
p = S5{r), (p) = {-iJjE,kJp\^E,kJ via 

(f'^ = a.i-cta.n{{py)/{p^)) + 7r6'(-(p^)) . (25) 

Due to translational invariance in a;-direction, (p) will only depend on the y-coordinate. 
For the following discussion, we will also discuss the at r = normalized expectation value 
((?) = {p)/{n) with (n) = {tpE,k^\S{r)\'ipE,k^) ■ This shall not be confused with the bulk 
expectation of a as it appears in the Hamiltonian. 

We will distinguish the two different cases of the half-plane y > (scattering from the 
lower or bottom boundary) and y < (scattering from the upper or top boundary). We 
shall further assume a plane wave with fc^, > moving in positive x-direction. The results 
for kx < are then obtained by changing bottom to top boundary and vice versa. The 
results for the K'-po'mt can also be deduced from the following discussion (see appendix 
A). The sign of A determines the sign of the expectation value of r and a. In the following, 
we set A = |A|, but in some of the following expression we explicitly use |A| for sake of 
clarity. 

We will discuss two different types of confinement. First, we use the fact that Dirac 
fermions can be confined by an infinite mass boundary, first discussed by Berry and Mon- 
dragon [57] . We then also study the reflection from a zigzag boundary first addressed in 
Ref. M- 



3.1 Infinite mass boundary 

With '4'E,k^ — ("01, "02, "03; "04)"^: the boundary conditions at the infinite mass boundary 
read (see appendix B and C) 



bottom 



lp4 



bottom 



■03 



top 



;02 

-04 



= -1 



(26) 



top 



Note that there are different boundary conditions depending on whether one approaches 
the boundary from below or above. 



3.1.1 Scattering behavior for plane waves of type I 

We first consider a plane wave scattered at y = with Ai — 1 and A2 = 0. The boundary 
conditions yield the following expressions for Ri, R2: 

2{ziCi±Si){Z2S2±C2) + {ziSi±Ci){Z2C2±S2)ziZ2 

= T^i? Z — T7 Z — ^ n Z — V Z — ^ — ^^'1 

(ziSi ± Ci)[Z2S2 ± C2)Zi + (ziCi ± Si){Z2C2 ± S2)Z2 
„ ^ 2 (ziCi ±Si)^ - (ziSi ±Ci)^zf 

-K2 — TZ21 T — v Z — ^ T? Z — V Z — ^ — y^^") 

(ziSl ± Ci){Z2S2 ± C2)Zi + (ziCi ± Si){Z2C2 ± 52)^2 

Above, we introduced the abbreviations = cos(i9^/2), — sin(z?^/2), and z^ = {kx + 
^ky)l \/^x + (^a)^i ^^ € {1,2}. The upper (lower) sign holds if the electron is scattered 
from the upper (lower) boundary. 

Let us first discuss the scattering behavior from the lower boundary. For kx = k cos Lp, 
the incident spin direction is given by Lpa = 7r/2 — \lp\. On the left hand side of Fig. [21 the 
reflected spin direction Lp'^ of Eq. (|25p is plotted against the incident spin direction Lp„. 
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reflection from botlom boundary 



rellection from top boundary 




incident spin direction (p^ [in n/2] incident spin direction (p^ [in n/2] 



Figure 2: The reflected versus the incident spin direction at y = for an incident plane wave 
with A2 = (type I) for various energies E. We use fivp = 5.6eVA and A = 6meV. Left: 
Reflection from the lower boundary. Right: Reflection from the upper boundary. 



At large energies with e — \/{hvpk) <C 1 and e <C sin^ (p, we have i?i = (1 — e) cosip 
and i?2 = iainLp — 2ecos(/9 and the spin polarization is approximately conserved. The 
expansion of Eq. ((25)) yields 

/ , cos Lp„ 

ip„^iPa+ <^-r—. • (29) 

1 + sm ip^ 

For energies close to the band gap energy of the type I-spinors, E 2A, scattering 
from the boundary acts as a spin polarizer since ip'^ ^ ipo ~ arctan(l/(2'y2)) ~ 19.5° for 
all incoming spin directions (pa- This angle corresponds to (ctj,) = 1/3. For E = A(2 + e^) 
with e <C 1, we get 

/ 1 \ 2 \/2 
ip'^ = arctan f ^ j + 3^ cos "/^ + (cos(2(^) - 5) . (30) 

This is a surprising result since i?i — > — 1 and incoming and reflected wave seem to compen- 
sate. But even though R2 — > — V6ee~*''' sin ip tends to zero, its admixture has a dominating 
effect. 

For the upper boundary, we obtain the expansion 

f I \ 2 , \/2 2 



ip'„^7T + arctan j - -e cos + — e (cos(2(^) - 5) . (31) 



Note that the different sign compared to Eq. ((30)) results in a different asymptotic behavior 
for large energies since f'^iE = 2A) is larger than the maximal incident spin direction 
ipa = TT. This different behavior is illustrated on the right hand side of Fig. O 

3.1.2 Scattering behavior for plane waves of type II with E >2X 

For a plane wave scattered at y = with Ai = Q and A2 = 1 with energy E > 2A, the 
boundary conditions yield the following expressions for i?i , R2: 

^ _p^2 (Z2S2±C2)^-(Z2C2±S2)^Z| 

^ (ziSl ± Ci)(z2S2 ± C2)zi + (ziCi ± Si){z2C2 ± 52)^2 

2{ziCi±Si){z2S2±C2) + {ziSi±Ci){z2C2±S2)ziZ2 , 

rL2 = T^2 7 1 w i ^ 7 1 w 1 N 

(ziSl ± Ci){Z2S2 ± 02)^1 + (ZICI ± Si){Z2C2 ± S2)Z2 
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Figure 3: The reflected versus the incident spin direction (rotated by tt) at ?/ = and = 
for various energies E > 2X. We use fivp = 5.6eVA and A = 6meV. Left: Reflection from the 
lower boundary. Right: Reflection from the upper boundary. 



For {E-2\\\)/{E + 2\\\) > {cosLpf, the abbreviations are the same as in Eqs. (I27p and 
(f28|) . For (£:-2|A|)/(i; + 2|A|) < (cosLpf, the reflected momentum kl = ±i q IS imaginary 
with 

fiVFq = ^-(-^-1-^1)^ + ^^ + {T^vfKY ■ (34) 

The sign is determined to yield an exponential decay in the reflected region. In Eqs. (I32p 
and ([33|) . zi is thus replaced by zi {kx T Q)/\/kx ^ Q^i where the upper (lower) sign 
holds for reflections from the upper (lower) boundary, and si by Si i-\/(cos'i?i — l)/2. 

Let us first discuss the scattering behavior from the lower boundary. On the left hand 
side of Fig. [3l the reflected spin direction is plotted against the incident spin direction 
rotated by tt. For large energies and normal incident ip « 7r/2, we again obtain ip'^ = Lpa- 
But for nearly parallel incident such that (E' — 2|A|)/(i? + 2|A|) < (cos</3)^, we obtain 
ip'^ — ±7r/2. For energies close to the band-gap E 2A, all reflected modes of type I are 
evanescent and scattering from the wall acts as a switch which leads to either ip'^ = ii/2 
or ip'^ = -7r/2. 

Let us understand the appearance of the two extreme values of ip'^ = ±7r/2 in the 
regime where ky is imaginary. Since zi is real and the incident and reflected wave of type 
1X2,+) compensate, the expectation value in a;-direction {a^) = 0. For the incident wave, 
( I CTy I) incident IS negative and for small incident angle, we thus have ip'^ = — 7r/2. But if |i?i| 
is large, the admixture of can lead to Lp'^ = n jl. Additionally, the spin in z-direction 

(az) assumes a non-zero value to guarantee |((t)| = 1. On the left hand side of Fig. [H this 
general behavior is shown whether the reflected spin angle (rotated by tt) , the expectation 
values ((Ti) (i = x,y,z) and the absolute value of the reflection amplitudes |i?i| and |i?2| 
is plotted versus the incident spin direction at y = at energy E = 3A. 

The scattering behavior from the upper boundary is considerably simpler. There, only 
two regimes appear with are marked by whether ky is real or imaginary. This can be seen 
on the right hand side of Figures [3] and ID 

3.1.3 Scattering behavior for plane waves of type II with E < 2\ 

For energies with E < 2A, one of the reflected modes becomes evanescent which leads to 
{(Jx) = 0. For a more detailed analysis, we have to distinguish the two cases E > X and 
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relleclion from bottom boundary rellection from top bou ndary 




incident spin angie (<Pg+7z) [in incident spin angle [in n/2} 

Figure 4: The reflected spin angle (rotated by tt), the expectation values (aj) {i = x,y,z) and 
the absolute value of the reflection amplitudes |i?i| and |i?2| versus the incident spin direction 
at y = and A2 = for energies E = 3A. We use hvp = 5.6eVA and A = 6meV. Left: 
Reflection from the lower boundary. Right: Reflection from the upper boundary. 



E < X. 

For X < E < 2X, the reflected momentum ky — ztiq is imaginary with the same 
expression as in Eq. ()34|) . The sign is determined to yield an exponential decay in the 
reflected region. With the ansatz 



+ i?ie-«l^l|xi.+ (fc.,±*g))+i?2e-^'='^|x2.+ (fc.,-fc^)) 



(35) 



we obtain the same expressions for Ri 
placement ci 



i?i and i?2 as in Eqs. ([5^ and ([55]) with the re- 



a/(1 + cosi?i)/2, si iA/(cosi9i - l)/2, and zi -i{ka: T q)/\/<l'^ - kl, 
where the upper (lower) sign holds for reflections from the upper (lower) boundary. 

Let us first discuss the lower boundary. For small incident spin direction, {ay) > and 
becomes zero at ipa- = fE < </?o ~ 19.5°. The reflected spin angle is thus Lp'^ — 7r/2 for 
if > (fiE and ip'^ = — 7r/2 for ip < LpE and for _E — A we have (fiE^x = 0. 

For the upper boundary, we have (ay) < for all angles and energies. In both cases, 
we have (ctz) ^ to fulfill the sum rule |((t)| = 1. 

For energies with < i? < A, there is no reflected wave of type I, 1x1,+), but one of 
the reflected momenta of |x2,+) is imaginary, ky — ztiq with the same definition as in Eq. 
With 



we have 



R2 

i?2 



\X2,+ {k.,±^q)) + R2e-<y\x2Ak., -k^)) 

{Z2S2 ± C2Y - (Z2C2 ± S2)^z| 

(Z2S2 ± £2) (2:202 ± S2)Z2 - {Z2C2 ± S2){Z2S2 ± 02)^2 
2 (i2S2 ± C2)(2:2S2 ± C2) - (z2C2 ± S2){z2C2 ± ^2)^2^2 
(Z2S2 ± C2)(2;2C2 ± 52)2:2 - (Z2C2 ± S2)(22S2 ± 02)^2 



Tzljrr-: 



Tz 



(36) 

(37) 
(38) 
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Figure 5: The reflected versus the incident spin direction at y = with A2 = (left hand 
side) and Ai = (rotated by vr) (right hand side) for various energies E >2X in the case of a 
zigzag boundary. We use Tivp = 5.6eVA and A = 6meV. 



with 52 = V(l + cos??^)/2, ~S2 = V(cos??i - l)/2, ~Z2 = -i{kx T q)l^/q^ - kl, and dl = 
|A|/(|A| — E). In the above equations, the upper (lower) sign holds for reflections from the 
upper (lower) boundary. 

We obtain {ay) = — 1 for the upper and {uy) = 1 for the lower boundary, respectively 
which is independent of the incident direction nor of the energy. 



3.2 Zigzag boundary 

Graphene can be terminated by a zigzag boundary which exposes only one sublattice to the 
boundary. With tpE.k^ — {ipi,ip2,'4'3,'4'4)'^ , the boundary conditions at a zigzag boundary 
thus read 

V'l = V'2 = (for bottom boundary) , tp^ = tp4 — (for top boundary) . (39) 

Here, we assumed that the bottom boundary is terminated by sublattice A and the top 
boundary by sublattice B. 

For a general plane wave Eq. ([23]) scattered at y ~ with energy E > 2X, the boundary 
conditions for the bottom boundary (sublattice A) Eq. ([39]) yield the following expressions 
for R2: 

2^l('Sl«2 + C1C2Z1Z2) + ^2(5202 - S2C2Z|) 

Hi = -Zi — (4Uj 

S1S2ZJ + C1C2ZIZ2 

2Ai{siCi - SiCizI) + A2{ciC2 + S1S2Z1Z2) 

H2 = ~Z2 (41 j 

C1C2Z2 + S1S2Z1Z2 

The boundary conditions for the upper boundary (sublattice B) yield the following ex- 
pressions for R2: 

2AiZi{CiC2 + S1S2Z1Z2) + A2Z2{S2C2- S2C2ZI) 

Kl = -Zi ■ (42) 

C1C2Z1 + SIS2Z2 

n _ 2AlZl{siCl- SlClzl) + A2Z2{siS2 + C1C2Z1Z2) , 

U2 — — Z2 : l^jj 

C1C2Z1 + S1S2Z2 
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incident of type I (A^=l A^=0) 



incident of type ti (A^=0, A^=1) 




200 400 600 800 1000 
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200 400 600 800 1000 
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Figure 6: Spin polarization in j;-direction as function of y for various energies E > 2\ with 
an infinite mass boundary. We use Tivp = 5.6eVA and A = 6meV. Left: Incident plane wave 
of type I. Right: Incident plane wave of type II. 



The abbreviations are the same as for the infinite mass boundary. Since the reflected angle 
is symmetric around normal incident, we will only discuss the reflection from the bottom 
boundary for kx > 

In Fig. [5l the reflected versus the incident spin direction at y = is shown for the two 
types of incident plane waves. As in the case of the inflnite mass boundary, {gx) = for 
incident plane waves of type II with cos^ ip > 2|A|)/(£' + 2|A|). But in contrary to the 
infinity mass boundary, the spin-polarization in out-of-plane direction assumes a non-zero 
value even when the reflected wave of type I is extended. For this case, i.e., ky S R, we 
obtain 

The A'-point yields the opposite sign such that there is no net-polarization in z-direction. 
For energies E < 2X a, similar discussion as in the case of inflnite mass boundary applies. 



4 Spin polarization close to the boundary 

So far we have only discussed the polarization properties at the boundary y = 0. For flnite 
y, we expect an oscillatory behavior of the reflected spin polarization. For E 2X and 
plane wave scattering of type I, fc^ ^ and the period will thus be solely determined by 
ky — > ^/2{2X/hvF). This oscillatory behavior is again independent of the incident spin 
polarization and results in a striped phase for the reflected spin-polarization. For E > 2A, 

1 /2 

two periods related to ky' contribute and a more complicated pattern emerges which also 
depends on the incident spin polarization and whether one deals with a reflection from 
the top or from the bottom. This hints to the fact that a Dirac particle in a box shows 
quasi-chaotic behavior [55]. 

In the following, we will study the spin polarization averaged over the incident direction 
for fixed Ai , A2 and including the two i^-points as function of the y-direction. We will 
further average over positive and negative /cjj-momenta. With an incident wave of type fj. 
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infinite mass boundary zigzag boundary 

n • I • I • I • I I • I • I • I • I r 




200 400 600 800 1000 200 400 600 800 1000 



y [in nm] y [in nm] 

Figure 7: Net spin polarization in x-direction as function of y for various energies E > 2A. 
We use hvF = 5.6eVA and A = 6meV. Left: Infinite mass boundary. Right: Zigzag boundary. 

and momentum k'^ = y/{E + (-1)^|A|)2 - X^/iTivp), [i e {1, 2}, we have 
1 ■ 1 

{p)^{r) = - - / d(p {ip E. ki- cos 5 {r~r)\-4)E.ki- cos v)k ■ (45) 

We only discuss the spin polarization at the lower boundary which depends on the sign of 
A (here we choose A = |A|). The spin polarization on the upper boundary is obtained by 
reversing the sign. 

In Fig. [6l the angle-averaged spin density A{px)'^{t') is shown as function of y for 
various energies E > 2X where A denotes the area of the sample. We show the results 
for an incident plane wave of type I (left hand side) and type II (right hand side) with an 
infinite mass boundary. There is a clear difference between the two types for low energies 
which is due to the appearance of imaginary momenta ky = zLiq for type Il-reflections. For 
low energies, most incident angles of the initial plane wave of type II lead to evanescent 
modes and thus to ((t^;) = 0. For large energies E > lO^X, the spin polarization of the two 
types have approximately the same absolute value, but differ in sign. 

Obviously, the above ensemble average breaks time-reversal symmetry since there is one 
incident plane wave with fixed fcj,-direction and two refiected plane waves. But if there is 
no coherence between the incident plane waves of type I and II, e.g., due to temperature, 
then time-reversal symmetry is effectively broken and we find a net polarization in x- 
direction by adding the two contributions (p)^ and (p)^^ (and possibly weighting them 
with the corresponding density of states). This is demonstrated in Fig. [71 where the sum 
of the two contributions AJ^fj.iP)'^ shown for a infinite mass boundary (left) and for a 
zigzag boundary (right). Moreover, we expect spin polarization in x-direction for various 
non-equilibrium situations. 

In the other two directions, we find no net spin polarization if the two inequivalent 
-ftT-points are included. We note, however, that py and pz assume a finite value for one 
-fC-point, only. This opens up the possibility of spin polarization in these directions in 
the presence of ripples or a magnetic field. Especially surface states due to, e.g., zigzag 
boundaries which effectively break the sublattice symmetry and which are not included in 
our continuous model should give rise to a finite spin polarization. 
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5 Dirac electrons with "pseudo-Rashba" spin-orbit cou- 
pling in nanoribbons 

In this section, we will consider graphene nanoribbons and the quantization properties of 
the transverse momenta in the presence of "pseudo-Rashba" spin-orbit coupling. We will 
then discuss the density and spin distribution at various energies. 



5.1 Quantization of the transverse momentum 

Let us first consider infinite mass boundaries. For a general plane wave with fixed momen- 
tum kx and energy E, ijjE,k^ {x, y) = (V'l, fp2, fp3, ^4)"^, there are four conditions that have 
to be satisfied, i.e., ipi = zLip^ and -02 = ±-04 at y = 0, and -01 = TV's ^nd V'2 = TV'4 at 
y = W, where the upper (lower) sign stands for the K{K')-pomt and W the width of the 
nanoribbon. For a zigzag nanoribbon which terminates on sublattice A at the bottom and 
on sublattice B at the top, the four conditions read = ■02 = at y = and ■03 = V'4 — 
a,ty = W. 

Let us first assume two propagating waves as in Eq. ((23)) . see also Fig. [51 In order to 
have a non-trivial solution, a necessary condition is 

det M = det ^ ^ ^ ) " {AB-'^ - AB'^) det B det B = , (46) 

with the bar denoting the complex conjugate. 

For infinite mass boundaries, the above matrices read at the ii'-point 

Si-CiZi C2-S2Z2 \ {si+ClZi)wi {C2+S2Z2)W2 



(ci-SiZi)zi {-S2+C2Z2)Z2 J ' V (ci + SiZi)ziWi -(52+0222)^:2^2 

(47) 

and for zigzag boundaries, we have 

.1 C2 \ ( C,Z,W, S2Z2W2 \ (^g^ 

CiZi —522^2 / \ SiZfWi ~C2Z2W2 J 

where we introduced — e^'^y^ and used the definitions of section [31 detM in Eq. 
is real and thus yields the quantization of the transverse momentum in y-direction. 
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For (hvpky)'^ < 4EX, there is the appearance of evanescent modes since ky = ztiq 
is imaginary. In this case, a general plane wave with fixed momentum and energy 
E > 2|A|, ipE.k.A^,y) = (V^i, V'2,V'3,V'4)^, is written as 

+ R^e"^y\xiAk.,iq))+R2e-''=ly\x2Mk.,-kl))\, (49) 

with 



hvpq = -{E - \\\)-^ + y + [fivk^Y , hvpkl = ^{E+\X\)^ - X'^ - {hvk.y . (50) 

Again, in order to have a non-trivial solution, Eq. (1461) must hold, but this time the 
matrices for infinite mass boundaries at the if-point read 

j^^f isi-ClZ^)wi C2-S2Z2 \ Sl+ClZ'^ {C2+S2Z2)W2 

\ {Cl - SlZi)z^Wl {-S2 + C2Z2)Z2 J ' \ {ci + SlZ^)z^ - (s2 + 0222)^2 W2 

(51) 

and for zigzag boundaries, we have 

A = ( ) ^B^( 7^,2 ~''T' ) , (52) 

\ CiZ{Wi ~S2Z2 ) \ Sx{zjY -C2Z^W2 J 

with wi = zf = {k^±q)/^kl-q\ ci ^ v/(l + cost?i)/2andsi -> i ^(cos i^i - l)/2. 

The definitions for the plane wave of type II remain unchanged. Since the wave function 
of the evancccnt mode is now real, the matrices A, B are not the complex conjugates of 
A, B, but given by 

Sl-CiZf C2 - S24 \ B-( (s1+CiZj7)wi {c2 + S2Z^)wl 



[ci - SlZ^)z^ {-S2 + C2Zl)zl J ' \^ [ci + SlZ^)z^Wl -{S2 + C2Zl){z2W2)* 

(53) 

for infinite mass boundaries, and for zigzag boundaries they read 

A^i ) , S ^ f -1 ) . (54) 

It is now preferable to write Eq. (|46p in powers of wi. For zigzag boundaries, this yields 
det M = 2ilm (dot B dot A) + 2wi{z'^ — z^){z2 — 22)siCiS2C2 (55) 



Zi {'W2Z2 - W2Z2)siCiS2C2 + (z^ ) {w2 - W2)S-^S2 + (W2-Z2 " W2Z2)c^C2 



0, 



which is purely imaginary and thus again yields a quantization of the transverse momentum 
in y-direction. For infinite mass boundaries, we obtain a similar expression. 



5.2 Spin and density distribution 

The particle density (n) at energy E is now obtained by summing over all transverse modes 
n that obey the above boundary conditions or the corresponding boundary conditions for 
the _ftr'-point. Denoting the rt-th transverse momentum of type II by A:^ „, we have 

(,i)(f) =Y^Y.ki^E,kMr- r')|V'£,feJ<5£.£,^,,.^ . (56) 

n 

In Fig. [9l the density distribution of a graphene nanoribbon of width W — lOOnm 
for various low energies with infinite mass (left) and zigzag (right) boundaries is shown. 
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infinite mass boundary zigzag boundary 

16 



14 

A 



0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

y/W y/W 

Figure 9: The density distribution of a graphene nanoribbon of width W = lOOnm for various 
low energies. We use hvp = 5.6eVA and A = 6meV. Left hand side: Infinite mass boundaries. 
Right hand side: Zigzag boundaries. The corresponding density distributions without "pseudo- 
Rashba" spin-orbit couphng are also shown (dashed lines). 




In general, the number of modes is the same with and without "pseudo-Rashba" spin- 
orbit coupling and the resulting density distributions only differ slightly. But for zigzag 
boundaries a,t E = 4A, we observe strong deviations due to the fact that there are 8 modes 
in the case with spin-orbit coupling in contrast to 12 modes in the case without spin-orbit 
coupling. Also note that whereas for the case without spin-orbit coupling, all modes are 
extended, some modes for the case with spin-orbit coupling are evanescent for the type 
I-branch. For zigzag boundaries, e.g., we have no extended and 4 evanescent (type-I) 
modes at E = 2X, 6 extended and 2 evanescent modes at E' = 4A and 6 extended and 6 
evanescent modes at E — 6A. 

The spin polarization at the boundaries is in all cases zero; in x-direction it is zero also 
for only one A'-point, in y- and z-direction it is non-zero for one A'-point, but averages to 
zero when two AT-points are included. This is an immediate consequence of time-reversal 
symmetry. In Ref. |30j . spin polarization in z-direction is reported for various kx-values 
within a lattice model of a zigzag nanoribbon. At equilibrium, this can only be attributed 
to edge-states which effectively break the sublattice symmetry and which are not included 
in our continuous model. 



6 Summary 

In this paper, we have investigated the spin dephasing of Dirac fermions with "pseudo- 
Rashba" spin-orbit coupling due to the reflection from a hard wall. In order to confine the 
Dirac electrons, we used infinite mass and zigzag boundaries. For large energies compared 
to the spin-orbit coupling, we obtained the expected result that there is hardly spin- 
dephasing due to the scattering process. But for energies close to the band gap for plane 
waves of type 1, E fa 2A, strong spin dephasing is observed. If the incident plane wave is of 
type II (gapless branch), even stronger effects are seen like the appearance of evanescent 
modes. We also observe the rotation of the spin in out-of-plane direction away from the 
boundary and for incident plane waves of type II also at the boundary. This effect will be 
canceled by averaging over the two inequivalent A'-points. 

We also discussed the spin polarization averaged over the incident direction and includ- 
ing the two A'-points. We find that for energies E > 2A, there is a finite spin polarization 
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in a;-direction when there is no coherence between the two branches. This polarization 
differs in sign for the upper and lower boundary, respectively. Also for non-equilibrium 
situations, there will be a spin polarization in this direction. 

We finally analyzed the spin and density distribution of graphene nanoribbons. At 
certain energies, the number of transverse modes does not match the one of a correspond- 
ing nanoribbon without "pseudo-Rashba" spin-orbit coupling. This results in significant 
changes in the density distribution. But generally, the "pseudo-Rashba" spin-orbit cou- 
pling leads to marginal differences, only. Further, there is no spin polarization if both 
K-points are included, but we find a finite spin polarization in y- and z-direction for one 
K-point, only. Surface states due to, e.g., zigzag boundaries which only live on one sublat- 
tice and thus break the valley-symmetry should therefore yield a finite spin polarization. 
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8 Appendix A: The full model including the two K- 
points 



The full model including the two iiT-points reads 

H = V {pxUzTx +PvTy) + A {KzTx<7y ~ Tyffx) , 



(57) 



where = ±1 denotes the two inequivalent if -points. For a given wave vector k the 
Hamiltonian around the if '-point {kz = — 1 ) reads 



n{k) 



/ 



Uvpikx — iky) 

y 2iX hvpikx — iky) 



hvF{kx + iky) —2iX \ 

hvpikx + iky) 



/ 



(58) 



The Hamiltonian around tlic^ if -point can thus be obtained from the Hamiltonian 
around the if -point by interchanging the pseudo-spin index and reversing the sign. All 
previous results without the mass term can thus be used. The results involving the mass 
term are obtained by M — » — M. This leads to a change in the boundary conditions, i.e., 



K' 



bottom 



3^2 
■^4 



K' 



bottom 



K' 



top 



1p4 



K' 



top 



(59) 



9 Appendix B: Massive Dirac fermions with "pseudo- 
Rashba" spin-orbit coupling 

Massive Dirac fermions with "pseudo-Rashba" spin-orbit interaction can be described by 

n = vp-T + X{TX a)-ez + Mv'^Tz, (60) 

where, among standard notation, A is the spin-orbit coupling parameter, and the Pauli 
matrices r, a describe the sublattice and the electron spin degree of freedom, respectively. 
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Squaring the Hamiltonian, we obtain the same eigenvectors as for massless Dirac 
fermions given in Eqs. (O and ©■ In the basis (|ai), |q;2), |/32)) the Hamiltonian 
reads 



with 
and 



n{k) = 



q± = ±hvF{kx ± iky)f±{\X\/hvFk) 



1 m 


1*+ 







q+ 


— m 














m 


q~ 


V 







—m J 



m — Mv 



f±{x) = Vl + ^±x. 
Again we find two types of solutions. The first type has eigenvalues 

ei,± = ±^M2i;4 + {hypk)^ + 2A2 + 2\\\^/(hvFky + 
with eigenspinors 



lxi,+(fc)> = 



/ sin(i?/2)cos(Ci/2) \ 
cos(i?/2)cos(Ci/2)e"' 
cos(\9/2)sin(Ci/2)e*'^ 

\ sin(i9/2)sin(Ci/2)e"'e*'/' / 



\xi.-m = 



(61) 

(62) 
(63) 

(64) 

\ 



with G [0,7r] and 

cosCi/2 = - , 

The second type has eigenvalues 



sin(t9/2)sin(Ci/2) 
cos(i?/2)sin(Ci/2)e"' 
-cos(z9/2)cos(Ci/2)e*'^ 
V -sin(z9/2)cos(Ci/2)e*''e*'^ / 

(65) 



_ fcz + iky 



e2,± = ± JM2u4 + {hvFky + 2A2 - 2\X\y/{hvFk)^ + X"^ 



(66) 



(67) 



with eigenspinors 

IX2,+ (fc)) 



IX2.-(fc)) 



/ cos(i?/2)sin(C2/2) \ 
-sin(z?/2) sin(C2/2)e*'' 
~sin(il/2)cos(C2/2)e»'^ 
V cos(T?/2)cos(C2/2)e*''e*''' / 

(68) 

Let us now consider expectation values within the eigenstates with wave functions 



/ cos(i?/2)cos(C2/2) \ 
- sin(??/2) cos(C2/2)e"' 
sin(i?/2)sin(C2/2)e*''' 
\ -cos(i?/2)sin(C2/2)e'''e*'^ / 



{r\k,fi,±) = --^\x^,±) 



/Lt e {1,2}, and A being the area of the system. Here we find 

sin I? sin cos cp 
{k, fi,±\f\k, fj,,±) =± I sin i9 sin sin iy9 



(69) 



(70) 



cosC 



and 



— sin i} sin ip \ 

(/c, l,±|CT|fc, 1,±) = -(A;,2,±|ct|/c,2,±) = ( siniJcosv? . (71) 

=Fcos?9cosCi/2 / 

Here we have assumed a positive spin-orbit coupling parameter, A = |A|, and is the usual 
azimuthal angle of the wave vector, k — k{cosip^ sin 1^9). Note that massive Dirac fermions 
assume a non-zero expectation value for the pseudo-spin and spin in z-direction. 
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10 Appendix C: Scattering from infinite mass bound- 
ary 

Dirac fermions can be confined by an infinite mass boundary, first discussed by Berry and 
Mondragon [27]. In the following, we will study the scattering behavior from a boundary 
located at y = and y = W . Within the strip < ?; < W , the mass of the Dirac fermions 
shall be zero; outside the strip, the mass shall be infinite. 

A general plane wave within the strip with fixed momentum fc^ and energy E > Q can 
be written as 

i^E,kAx.y) = e^''-^\A,e<y\xi,+ {kx,kl)) + A2e<y\x2Ak..kl)) 
+ R^e-'Ky\xiAk.,-kl)) + R2e-<y\x2Ak.,-kl)) 



with 



nvpk^; = ^{E + {-iY\\\f-x^-{hvFk^Y 



(72) 



(73) 



Aie{i,2}. 

The wave function of the transmitted electron is also decomposed by the two eigen- 
functions 



i^E,kAx.y) = e''^-" \T^e<y\xiAk., kl)) + T2e''=2^|x2,+ (fc., fc^)) 



with 



nvpk'^ ^ \J{\/e'^ - Aih^+{-iY\\\Y - A2 - {hvpk^y 

In the limit M oo, the transmitted plane wave simplifies 



1 



e ±2 



( 1 \ 

1 

V / 



(74) 
(75) 

(76) 



and 



o 1 1 



e ±2 



( 1 \ 

-1 

V -^A J 



(77) 



/ 1 \ 

-1 
V ^A J 

with sx — A/|A|. The different expressions at y = and y — W originate from the 
different sign of fivpky — > zLiMv"^ that has to be chosen to yield an exponential decay in 
the infinite mass region. It therefor only depends on whether one deals with the upper or 
lower boundary. 

At the boundaries y — and y — W, the four components have to be continuous to 
guarantee a continuous current which leads to the following two sets of equations: 



■^E.k, [x, W) = IpE.k^ [x, W) 



(78) 



With il)E,k^ — (V'lj "02, V'Sj i^iY' ^ this translates to the familiar boundary condition from 
Ref. [57] for the two spin channels, respectively: 



y=0 



V'4 



■03 



y=W 



;02 

04 



(79) 



y=W 
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